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ABSTRACT 

We have discovered decimetric signatures of the chromospheric evaporation process. Evidence for the radio 
detection of chromospheric evaporation is based on the radio-inferred values of (1) the electron density, (2) the 
propagation speed, and (3) the timing, which are found to be in good agreement with statistical values inferred 
from the blueshifted Ca xrx soft X-ray line. The physical basis of our model is that free-free absorption of 
plasma emission is strongly modified by the steep density gradient and the large temperature increase in the 
upflowing flare plasma. The steplike density increase at the chromospheric evaporation front causes a local 
discontinuity in the plasma frequency, manifested as almost infinite drift rate in decimetric type III bursts. The 
large temperature increase of the upflowing plasma considerably reduces the local free-free opacity (due to the 
T~ 3/2 dependence) and thus enhances the brightness of radio bursts emitted at the local plasma frequency 
near the chromospheric evaporation front, while a high-frequency cutoff is expected in the high-density regions 
behind the front, which can be used to infer the velocity of the upflowing plasma. From model calculations we 
find strong evidence that decimetric bursts with a slowly drifting high-frequency cutoff are produced by funda- 
mental plasma emission, contrary to the widespread belief that decimetric bursts are preferentially emitted at 
the harmonic plasma level. 

We analyzed 21 flare episodes from 1991-1993 for which broadband (100-3000 MHz) radio dynamic 
spectra from Phoenix, hard X-ray data from BXTSE/CGRO, and soft X-ray data from GOES were available. 
We detected slowly drifting high-frequency cutoffs between 1.1 and 3.0 GHz, with drift rates of —41 + 32 
MHz s" 1 , extending over time intervals of 24 + 23 s. Developing a density model for type Ill-emitting flare 
loops based on the statistically observed drift rate of type III bursts by Alvarez & Haddock (1973), we infer 
velocities of up to 360 km s -1 for the upflowing plasma, with an average of u CB = 236 ± 130 km s _1 for 
episodes with 5-15 s duration. The mean electron density of the upflowing plasma is n e = 5.2( + 3.1) x lO 10 * 
cm -3 when it is first detected in radio, at coronal altitudes of h 0 = 9.2 + 2.3 Mm. 

Subject headings: Sun: chromosphere — Sun: flares — Sun: radio radiation 


1. INTRODUCTION 

In the early phase of a solar flare, the plasma in a flare loop 
displays dynamic processes such as turbulent motion (with 
velocities exceeding 100 km s -1 ) and high-speed plasma 
upflows (with bulk velocities of 300-400 km s _1 ), as inferred 
from the line broadening and blueshift of soft X-ray (SXR) lines 
in Ca xix and Fe xxv (Antonucci et al. 1982; Antonucci, 
Gabriel, & Dennis 1984). The upward motion of the flare 
plasma in the early impulsive phase is generally referred to as 
“chromospheric evaporation” (Sturrock 1973) and is believed 
to be a consequence of local heating of the chromosphere near 
the footpoints of flare loops, produced either by collisions from 
precipitating electrons (Canfield et al. 1980) or by heat conduc- 
tion (Antiochos & Sturrock 1978). Chromospheric evaporation 
is considered to be the main mechanism for transporting the 
hot SXR-emitting flare plasma to coronal levels (Antonucci et 
al. 1984). Reviews on this subject can be found in Doschek et 
al. (1986), Antonucci (1989), Doschek (1990), and Antonucci et 
al. (i" 4 ). 

The chromospheric evaporation process has thus far been 

studied chiefly in SXRs and Ha. In this paper we address for 
the first time observational evidence at radio wavelengths. The 


basic idea is that the local disturbance of the electron density 
and temperature, introduced by the upflowing chromospheric 
plasma, is detectable from radio bursts emitted at the local 
plasma frequency. Plasma emission produced by electron 
beams has been observed in the lower corona up to a frequency 
of 8.4 GHz (Benz et al. 1992). The detection of plasma emission 
at such high frequencies requires overdense flux tubes, so that 
plasma emission can escape in a direction perpendicular to the 
flux-tube axis, where the density scale height is much shorter 
than in a homogeneous corona, and thus, free-free absorption 
is substantially reduced. In the event of chromospheric evapo- 
ration, we expect that the upflowing plasma surrounds over- 
dense flux tubes and seals off - escape routes for plasma 
emission, because the additional plasma material, if sufficiently 
dense, makes the escape routes optically thick owing to free- 
free absorption. Since the evaporating plasma propagates 
upward with a bulk speed of «300 km s~\ it is expected to 
produce a slowly drifting high-frequency cutoff for plasma 
emission. This high-frequency cutoff is thought to apply to any 
kind of plasma emission originating in “evaporating” flare 
loops, e.g., to type III bursts excited by precipitating electron 
beams. The drift rate of this high-frequency cutoff for plasma 
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TABLE 1 


Parameters of Plasma During Chromospheric Evaporation 



Velocity 

Temperature 

Density 

Plasma Frequency 


V CE 

T e 

(MK) 

”e 

v p 

Time Phase 

(km s" 1 ) 

10 11 cm" 3 

[GHz] 

Flare start ...... 



0.56-0.96 

2.1-2.8 

Flare end 



2.0-3.5 

4.0-5.3 

Average 

270 ± 90 

16.5 ± 2.4 

0.8-1.5 

2.5— 3.5 


emission is estimated to be 


/ dv\ dv dh 

\dt/ CE dh 8t 



cos S , 


( 1 ) 


where v is the observed frequency (assumed to be at the funda- 
mental plasma frequency), X the local density scale height, h the 
altitude, p CE the velocity of the chromospheric evaporation 
front, and # the propagation angle with respect to the vertical. 
Typical physical parameters of the chromospheric evaporating 
plasma are listed in Table 1 (from Antonucci et al. 1984). Based 
on these electron densities, we calculate that the plasma fre- 
quency of the upflowing plasma varies in the range of v p = 
2.5— 3.5 GHz and can be as low as 2.1 GHz at the start of the 
flare. For an upward moving (5 = 0) evaporation front, a 
velocity of v CE « 270 kms _1 ,a scale height range of X = 5-66 
Mm (inferred from decimetric type III bursts; Aschwanden & 
Benz 1986), and a frequency of v = 2 GHz, we estimate a drift 
rate of ( dv/dt) CE = —8 to —108 MHz s" 1 for the high- 
frequency cutoff. From this consideration we expect slowly 
drifting high-frequency cutoffs at decimetric frequencies for 
various radio bursts that are related to plasma emission, pref- 
erentially at the start of the impulsive phase of flares. 

This paper is dedicated to testing these theoretical predic- 
tions and to exploring the effects of chromospheric evapo- 
ration at decimetric wavelengths. We search for evidence of 
radio signatures produced by chromospheric evaporation by 
using decimetric observations from Phoenix, which has been 
recently equipped with broadband receivers in the extended 
frequency range of 1-3 GHz. A representative selection of 
radio bursts observed in this little-explored frequency range 
can be found in Isliker & Benz (1994). A substantial number of 
those radio bursts show slowly drifting high-frequency cutoffs 
at frequencies > 1 GHz, as expected from the considerations 
above. 

In § 2 we develop some theoretical tools to analyze deci- 
metric signatures of chromospheric evaporation, in § 3 we 
analyze flare observations at decimetric wavelengths and in 
X-rays, in § 4 we discuss the results, and in § 5 we summarize 
the findings and conclusions. 


2.1. Flare Scenario 

A general scenario for the chromospheric evaporation 
process is depicted in Figure 1, where the typical time evolu- 
tion of radio, hard X-ray (HXR), and SXR emission during 
different flare phases are outlined. Assuming that the primary 
energy release site where particle acceleration occurs is located 
in the corona, electron beams can be observed traveling both 
upward (as metric type III bursts) and downward (as deci- 
metric reverse-slope [RS] drifting bursts) during the impulsive 
flare phase and sometimes in the preflare phase before HXR 
emission is detectable (Fig. 1, top). At about the onset of HXR 
emission, a chromospheric evaporation front starts to rise from 
the footpoints of the flare loops, characterized by a high tem- 
perature and a steep density gradient at the front. While the 
evaporation front ascends in the flare loops, it will intercept 
with downward-propagating electron beams (denoted RS in 
Fig. 1) in the flare loops. At the interface between the ascending 
evaporation front and the cooler flare loop, the free-free 
opacity will be abruptly reduced due to the higher temperature 
of the upflowing plasma, while the opacity increases again 
behind the evaporation front due to the increasing density. 
Thus, decimetric type III bursts intercepting the evaporation 
front will be observed at a higher brightness temperature at the 
interface of the evaporation front due to the reduced free-free 
opacity, and will exhibit a much higher frequency-drift rate due 
to the locally steep density gradient at the evaporation front. 
After the (RS) type III bursts have crossed the evaporation 
front, plasma emission will be completely absorbed due to the 
higher density and increased path length in the high-density 
medium, which is observed as a high-frequency cutoff. This 
cutoff drifts slowly toward lower frequencies (Fig. 1, top ) as the 
evaporation front reaches higher altitudes (Fig. 1, bottom ), 
because the local plasma frequency decreases with height. SXR 
emission is produced in the volume of the upflowing hot 
plasma, leading to a monotonic increase of the SXR emission 
measure during the impulsive phase (Fig. 1, middle). The 
trapped particles in the flare loop undergo kinetic plasma 
instabilities, occasionally resulting in coherent radio emission 
(e.g., observed in the form of decimetric radio pulsations). If the 
coherent radio emission is produced by waves that are related 
to the plasma frequency, e.g., harmonic plasma waves or upper 
hybrid waves, we expect that they are subject to the same 
high-frequency cutoff as the decimetric type III bursts (radio 
emission from the trap is marked with a hatched area in Fig. 1, 
top). The boundaries of the trap at the mirrorpoints are then 
pushed together from the ascending chromospheric plasma 
(Fig. 1, bottom). Eventually, the evaporated chromospheric 
material reaches the apex of the flare loops or reaches pressure 
equilibrium, and the high-frequency cutoff for plasma emission 
will disappear. 


2. THEORY 

First we need to develop some theoretical tools to perform a 
quantitative analysis of the data and to establish evidence for a 
high-frequency cutoff of plasma emission produced by the 
chromospheric evaporation process. In the following we 
outline a relevant flare scenario (§ 2.1), we derive a suitable 
density model for flare loops (§ 2.2), we determine the velocity 
of chromospheric evaporation fronts (§ 2.3), and perform 
numerical calculations of the free-free absorption for plasma 
emission originating in the environment of a chromospheric 
evaporating plasma (§ 2.4). 


2.2. Density Model for T ype Ill-Emitting Flare Loops 
The coronal electron density and temperature in flaring 
active regions are highly inhomogeneous. Thus, there is no 
standard density model that characterizes flare loops and the 
quiet corona equally well. Existing density models of the 
corona are based on scattered white light for the quiet corona 
(e.g., Baumbach-Allen model, Newkirk model), with an arbi- 
trary adjustment factor of 5-10 for active regions. The models 
assume hydrostatic equilibrium, either for an isothermal atmo- 
sphere (barometric model) or including the temperature gra- 
dient in the upper transition region by constraints of the 
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Preflare Impulsive Phase Postflare 


Fig. 1. — Flare scenario illustrating the evolution of chromospheric evapo- 
ration (CEF), nonthermal hard X-ray emission (HXR NT ), thermal hard X-ray 
emission (HXR Th ), soft X-ray emission (SXR Th ), and radio emission (type III 
bursts, reverse-slope [RS] bursts, and plasma emission from a trap). ACC 
denotes the coronal acceleration region. In the radio dynamic spectrum (top) a 
slowly drifting plasma frequency cutoff is expected owing to the density 
enhancement conveyed by the chromospheric evaporation front; CEFI and 
CEF2 mark the plasma frequencies corresponding to the low-density (CEFI) 
and high-density edge (CEF2) at the evaporation front. The density of the 
upflowing chromospheric plasma is indicated with contour lines (bottom), 
marking a steep gradient at the evaporation front that corresponds to the jump 
of the plasma frequency from CEFI to CEF2 (top). This regime CEF1-CEF2 is 
also the frequency range where the opacity changes for plasma emission are 
most dramatic, and where the drift rate of type III bursts becomes almost 
infinite. 


conductive flux (Dulk et al. 1977). The Baumbach-Allen model 
and the barometric model for an isothermal corona with a 
temperature of 1.5 MK are shown in Figure 2. These models 
are not appropriate for type Ill-emitting flare loops for two 
reasons: (1) they represent lower limits on the electron density 
in type Ill-emitting flare loops, because the escape condition 
for plasma emission (to overcome free-free absorption) dictates 
overdense loops (e.g., by a factor of % 18 in Benz et al. 1992), 
and (2) flare loops are subject to substantial heating and mass 
flow, and are therefore not in hydrodynamic equilibrium. 

Here we are interested in estimating the electron density in 
flare loops where chromospheric evaporation and plasma 
emission from injected electron beams occurs. We construct an 
empirical density model n e (h) by using the statistically observed 
drift rate of type III bursts according to Alvarez & Haddock 
(1973), 


dv 

dt 


- Av a , 


( 2 ) 


which follows a power-law function of the frequency, A and a 
being constants. At large coronal heights, hydrostatic equi- 
librium is considered to be a good approximation and a baro- 
metric density model can be justified with a density scale height 
corresponding to the coronal temperature. Thus, we combine 
the two models by assuming a power-law function in the lower 



Fig. 2. — Density model (thick line) for type Ill-emitting (overdense) flux 
tubes, consisting of a barometric component for h > h lt a power-law com- 
ponent for h < h u and a (time-dependent) enhanced density component 
associated with the upflowing chromospheric plasma (indicated here at a 
height of <;10 Mm). The Baumbach-Allen model (dashed line) matches the 
barometric model (thin line) for heights h < h v The plasma frequency scale is 
indicated on the right. 


corona and an exponential form in the upper corona, i.e., 


n e {h) = 


n l (h/h l y p 

n Q exp ( — hjX) 


h < h l , 
h>h 1 , 


( 3 ) 


with a transition at height h u specified by the requirement of a 
smooth transition between the two regimes. The derivation of 
the numerical constants n u h lt n Q and A is given in Appendix A. 
This density model for flare loops is shown in Figure 2. Since 
decimetric type III bursts occur preferentially in closed (flare) 
loops (Aschwanden, Benz, & Montello 1994), while metric type 
III bursts generally originate from electron beams propagating 
along open field lines, our density model represents a gradual 
transition from closed to open magnetic structures with height. 
Nevertheless, it is intended to provide a representative density 
function for type Ill-emitting structures in a statistical sense. It 
provides the basis for measuring speeds of the chromospheric 
evaporation front with radio methods. 


2.3. The Velocity of the Chromospheric Evaporation Front 

The drift rate of a type III burst is proportional to the exciter 
speed for a given frequency v (see eq. [A6] in Appendix A). This 
proportionality is equally true for any propagating agent that 
affects plasma emission. In particular, we can relate the slowly 
drifting high-frequency cutoff of plasma emission to the veloc- 
ity of the chromospheric evaporation front in the same way. 
Since both the type III bursts and the plasma frequency cutoff 
occur in the same flare loop described by our density model 
(eq. [3]), we can apply the drift rate relation (A6) to both 
exciters. Hence, the two relations ( dv/8t) m oc v m and 
{dvjdt) CE oc v CE have the same proportionality constant (in eq. 
[A6]), and the velocity r CE of a chromospheric evaporation 
front can be found from 


V CE 


= ”iij 


(dv/dt) CE 

(d\/dt) m ' 


( 4 ) 


The drift rate of type III bursts, (dv/dt) can be inserted from 
equation (2) of Alvarez & Haddock (1973). A representative 
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TABLE 2 

Type III Exciter Speeds 


Frequency Number of 

Velocity Range Analyzed 

v m jc (MHz) Events Reference 


0.45 45-60 MHz 8 1 

0.33 7 220 MHz 50 2 

0.15 ±0.03 160-320 MHz 12 3 

0.16+0.05 56 kHz 3 4 

0.14 + 0.04 1 MHz 28 5 

0.18 + 0.04 1.446 GHz 1 6 

0.15 ± 0.02 250-850 MHz 1 7 

0.30 ±0.10 560-1240 MHz 1 7 


0.15 ±0.04 56 kHz-1.5 GHz ... 8 


References. — (1) Wild, Sheridan, & Neylan 1959; (2) Stewart 1965; 
(3) Fokker 1970; (4) Lin et al. 1981, 1986; (5) Dulk et al. 1987; (6) 
Aschwanden et al. 1992; (7) Aschwanden 1993; (8) Weighted average 
(References after 1970). 


value of the type III exciter speed is v m jc = 0.15 ± 0.04 (see 
average value in Table 2, a discussion in the context of our 
model is given in Appendix B). Figure 3 shows the drifting 
plasma frequency cutoff in a radio dynamic spectrum, calcu- 
lated for chromospheric evaporation fronts with upward veloc- 
ities of v CE — 100, 200, and 300 km s“ l . Note that the detection 
of a plasma frequency cutoff at <3 GHz is delayed by more 
than 20 s with respect to the star of a chromospheric evapo- 
ration, for propagation speeds of v CE < 300 km. 

2.4. The High-Frequency Cutoff of Plasma Emission 
We calculate the radio spectrum of plasma emission in the 
environment of a chromospheric evaporation front, with par- 
ticular emphasis on the high-frequency cutoff. The free-free 
absorption coefficient in typical coronal conditions can be 
approximated (Dulk 1985) by 

ic ff «0.2^T- 3/2 v~ 2 cm” 1 . (5) 

In the case of plasma emission the observed frequency v can be 
expressed as a function of the local density n e by the plasma 



Fig. 3. — Evolution of the plasma frequency cutoff calculated for an 
upflowing chromospheric evaporation front with a constant velocity of v CE = 
100, 200, or 300 km s" \ according to the density model illustrated in Fig. 2. 


frequency relation (A3). The resulting optical depth along the 
path r is then 


Tf f (r) = 


K {( dr — 2.5 x 10 


-9 


n e (r)T(r)~ 312 dr . (6) 


The observed radio brightness is reduced by the absorption 
factor e~ xu . If the intrinsic brightness temperature of the radio 
source is characterized by the spectrum T B0 (v), the observed 
spectrum will be 


TO = TO) exp-^ lvl) , (7) 


where the frequency dependence of the location r(v) is given by 
the plasma frequency relation v(r) = s 8980 n e {r)~ 112 and the 
density model n e (r). 

In order to study the effects of a chromospheric evaporation 
front on the observed radio spectrum of plasma emission, we 
perform numerical calculations of the radio brightness tem- 
perature T fl (v). We consider the situation of a vertical loop 
embedded in the ascending chromospheric evaporation 
plasma. The electron density in the overdense loop is specified 
by our density model (eq. [3]) in the vertical direction, and by a 
Gaussian distribution with half-width d in the horizontal direc- 
tion. The chromospheric evaporation plasma is approximated 
by a horizontally stratified structure, which has a density that 
is enhanced by a factor of q CE = 4 relative to the embedded flux 
tube (suggested by the average density variation measured 
from Ca xix during a flare, see Table 1). The front of the 
upflowing plasma moves with a velocity of y CE and the inter- 
face is located at height h CE = v CE t for times t > 0. The two- 
dimensional distribution of the electron density, 


n Q exp ( - - ) exp 


-(x/d ) 2 


n e {h, x, t) = < 


e h 
n n \ -t 


p X 


exp 


(x/d ) 2 


h > /ij , 

h CE (t) <h<h lf (8) 
h < ^ce(^) > 


is shown in the form of a contourplot in Figure 4 (left). The 
structure of the magnetic field that defines the vertical flux tube 
is assumed to be unaffected by the upflowing plasma. For the 
temperature distribution we simply assume a coronal tem- 
perature of T c = 1.5 MK for the overdense loop and a hotter 
temperature of 7^ E = 15 MK for the chromospheric evapo- 
ration front, 


T(h, x, t) — 



h > h CE (t ) , 
h < ^ce(0 • 


(9) 


For a given line-of-sight with aspect angle 9 to the vertical, the 
coordinates are x = r sin (5) and y — h + r cos (9), for a ray 
starting at height h at the center (x = 0) of the cylindrical flux 
tube. The radio brightness spectrum T b [v(/j)] originating from 
the center of the overdense flux tube can then be calculated in a 
straightforward manner by numerical integration along a 
selected line of sight using the two-dimensional distributions 
n e (h, x, t) and T(h, x, t) (eqs. [8] and [9]) for a given time t. 

The resulting radio spectrum T B (v)/T B0 (v) is shown in Figure 
4 ( bottom right ) for a flux-tube radius of d = 2 Mm, for a height 
of h CE = 20 Mm of the chromospheric evaporation front, and 
for an aspect angle of >9 = 45°. Note that the spectrum is nor- 
malized by the intrinsic brightness temperature T B0 that corre- 
sponds to the spectrum without the effect of free-free 
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Plasma frequency [GHz] 



Fig. 4. — Two-dimensional density model with a vertical flux tube (with FWHM of 2 Mm) and a horizontally stratified chromospheric evaporation front (CEF) at 
a height of < 20 Mm. The two-dimensional density distribution is shown (left) with logarithmic contour levels (in steps of 0.25). The vertical density distribution is 
identical with the density model in Fig. 2, except that here the height of the CEF is 20 Mm. The dependence of the plasma frequency vs. height is shown ( top right). 
The radio brightness spectrum T B (v)/T B0 (v) of fundamental (s = 1) and harmonic (s = 2) plasma emission excited by an electron beam propagating vertically along the 
overdense flux tube is calculated for an aspect angle of 45° (bottom right). 


absorption (as defined in eq. [7]). The effect of the chromo- 
spheric evaporation front on the radio spectrum can clearly be 
seen in this example: (1) since the density structure at the 
evaporation front was assumed to be a step function with an 
enhancement factor of q CE = 4, the plasma frequency jumps by 
a factor of 2 (from v CE1 = 0.7 GHz to v CE2 = 1.4 GHz), while 
the opacity is almost constant over this frequency range; (2) 
above the frequency v CE2 , the brightness first increases, with a 
peak at «1.7 GHz, because the opacity (proportional to 
n e T~ 3/2 ) is strongly reduced due to the higher temperature 
(T ce /T c = 10); (3) above v » 1.7 GHz the radio spectrum 
rapidly falls off because the opacity gradually increases due to 
the increasing path length dr(h ) = ( h CE — h)/cos (9) in the high- 
density region behind the chromospheric evaporation front 
(h < h CE ). The transit of the chromospheric evaporation front 
produces a relatively narrowband structure in the decimetric 
spectrum, with a bandwidth of Av « 0.6 GHz at v » 1.7 GHz 
(i.e., Av/v « 0.35 in this example). For harmonic plasma emis- 
sion (dashed spectrum in Figure 4, bottom right), the effect of 
the hot evaporation front is much less dramatic and does not 
produce a strong change in the 1-3 GHz frequency range of 
interest here. The high-frequency cutoff is relatively broad and 
extends to > 10 GHz. Flux tubes with very large diameters 
(;>50 Mm) would be required in order to produce a similar 
effect in the < 3 GHz range for harmonic plasma emission. 

Time-dependent simulations for chromospheric evaporation 
fronts (i.e., h CE = v CE t) with different temperatures (T CE = 10 
MK and T CE = 20 MK) are shown in Figure 5, assuming a 
constant velocity of the evaporating plasma (u CE = 300 km 
s~‘), illustrated for a periodic sequence of bursts that emit at 


the fundamental plasma frequency. The chromospheric evapo- 
ration front passes the 3 GHz plasma level between 20 s and 36 
s and makes the flux tube optically transparent after this time, 
because of the higher temperature of the evaporating material. 
A new decimetric component becomes visible after 36 s, whose 
peak and start frequencies drift slowly toward lower fre- 
quencies. For the 20 MK hot chromospheric evaporation front 
the peak frequency is 2 GHz about 60 s after start, and drifts 
toward 1 GHz after 110 s. The start frequency of this deci- 
metric component drifts similarly toward lower frequencies, 
while the spectral cutoff at the high-frequency side is steeper for 
the cooler {T CE — 10 MK) chromospheric material. According 
to the theoretical calculations above, the variation of the radio 
start frequency is merely a temperature effect: the free-free 
opacity is drastically reduced after the arrival of the hot 
upflowing plasma due to the T _3/2 dependence in the free-free 
absorption coefficient (eq. [5]). Thus, the start frequency of 
radio bursts emitted at the plasma frequency bear an impor- 
tant potential for temperature diagnostics. 

We investigate the dependence of the radio spectrum on 
various parameters in Figure 6. The calculated radio spectrum 
T B (v)/T B0 (v) is shown for fundamental (s = 1 ; Fig. 6, left side ) 
and harmonic (s = 2; Fig. 6, right side ) plasma emission, for 
several values of the temperature of the upflowing plasma 
(7^ E = 7.5, 15, 30, MK; Fig. 6, top row), the aspect angle 
(9 - 10°, 45°, 80°; Fig. 6, middle row), and the flux-tube radius 
(d = 1, 2, 4 Mm; Fig. 6, bottom row). From these computations 
we see that (1) the fundamental plasma emission has a steep 
high-frequency cutoff which is below 3 GHz in most of the 
cases, while the harmonic plasma emission cutoff is more 
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Fig. 5. — Simulated evolution of the radio dynamic spectrum during the first 3 minutes of a flare using the scenario described in the text. The radio spectrum is 
calculated for a periodic sequence of bursts emitting at the fundamental plasma frequency and propagating along the center of an overdense loop of radius 2 Mm. 
The curved trajectories mark the plasma frequency at the low- and high-density edges of a chromospheric evaporation front propagating with a velocity of 300 km 
s _1 upward in the corona, with a temperature of 20 MK (top), or 10 MK (bottom). The density enhancement in the evaporation front is assumed to be a factor of 4. 
Note the enhanced radio brightness and the high-frequency cutoff associated with the (slowly drifting) passage of the evaporation front. 


smooth and is above ss 3 GHz; (2) a higher temperature of the 
upflowing plasma shifts the start frequency to higher values; (3) 
radio bursts at the solar limb (5 > 80°) have a steeper spectral 
cutoff and a lower cutoff frequency than do bursts from inter- 
mediate (3 as 45°) longitudes; (4) radio bursts from near the 
center of the solar disk are also subject to strong absorption 
because the column depth of free-free absorption along the 
vertical axis of the flux tube is large; (5) the flux-tube diameter 
essentially determines the suppression of plasma emission at 
the low-frequency side, while the high-frequency cutoff is not 
affected by the flux-tube diameter. 

3. OBSERVATIONS 

We analyze radio observations from the radio spectrometer 
Phoenix at Bleien, operated by ETH Zurich (Switzerland), 
HXR observations from the Compton Gamma-Ray Observatory 
( CGRO ), and soft X-ray observations from the Geostationary 
Operational Environmental Satellites (GOES). 

3.1. Instruments and Data 

The frequency-agile broadband spectrometer Phoenix, 
which covers the frequency range of 0. 1-3.0 GHz, was put into 
operation in 1989 June. Instrumental descriptions and first 


results are given in Benz et al. (1991) and Isliker & Benz (1994). 
Since 1992 September, Phoenix has been operated mostly in a 
broadband mode with continuous coverage over the full fre- 
quency range of 0. 1-3.0 GHz, using a time resolution of 0.1 s. 
These broadband observations are the most useful ones for the 
purpose of our study. We analyzed 30 events with simulta- 
neous detection of HXR emission by CGRO, plus an additional 
flare without CGRO data. In at least 11 out of these 31 flares 
we found radio signatures of chromospheric evaporation, 
which is a lower limit, since not all events had optimum fre- 
quency and time coverage in radio. 

Hard X-ray data from CGRO are obtained from the Burst 
and Transient Source Experiment (BATSE). We use BATSE 
background data (DISCLA data type), with 1.024 s time 
resolution, from the eight large-area detectors (LADs), with 
2000 cm -3 detector area each, in the four energy channels of 
25-50 keV, 50-100 keV, 100-300 keV, and >300 keV. Instru- 
mental descriptions of BATSE data can be found in Fishman 
et al. (1989, 1992). The BATSE data were provided by the Solar 
Data Analysis Center (SDAC) at NASA/Goddard Space Flight 
Center (GSFC). 

SXR data were obtained from the NOAA-operated GOES 
satellites, with wavelength ranges of 0.5-4 A and 1-8 A, corre- 
sponding to energy bands of 3-25 keV and 1.5-12 keV. The 
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Fundamental (s=1) plasma emission Harmonic (s=2) plasma emission 






Fig. 6. — Calculated radio spectra in parameter space, for fundamental (s = 1, left) and harmonic ( s = 2, right) plasma emission, and for different temperatures (T e , 
top row), aspect angles (9, middle row), and flux tube diameters (d, bottom row). The unvaried parameters are the same as in Fig. 4:d = 2 Mm, Ji CE = 20 Mm, T ce =15 
MK, q EV = 4, 9 = 45°, and a = 1.84. Note that the frequency scales are linear but different for fundamental (0.1-3.0 GHz) and harmonic (0.1-10 GHz) plasma 
emission. 


time resolution of the GOES data is 3 s. The SXR data were 
provided through the GOES Archive at SDAC. 

3.2. Data Analysis 

In the analysis of the radio data, we focus on those flare 
episodes in which a high-frequency cutoff is apparent in the 
radio dynamic spectra of 0. 1-3.0 GHz. We select time intervals 
with a continuous high-frequency cutoff", regardless of whether 
the enhanced emission below the cutoff frequency consists of 
continuum emission (e.g., type IV bursts) or of fast-drift bursts 
(e.g., type Illdm bursts or decimetric broadband pulsations). 
The analysis presented here contains 21 episodes extracted 
from 11 different flares. 

During a selected time interval t h i = 1, ..., n with n data- 
points (measured with a typical time resolution of 0.1 s) we 
determine the cutoff frequencies v h i — 1, ..., n (or start 
frequency) from the highest frequency where the radio spec- 
trum E/v, t = t t ) exceeds a threshold of 3 times the rms noise 
level above the background. In the case of bursty emission, we 
sample the median values of the start frequency v ; from larger 
time bins (with « 1 s duration). Almost all of the episodes show 
a high-frequency cutoff v ,{(,-) that drifts toward lower fre- 
quencies, while drifts toward higher frequencies were found to 
be rare, and are omitted in the subsequent analysis. 

We fit a theoretical dynamic spectrum v(t) of an exciter 
moving with constant velocity v toward a lower density region. 
Using the relation for the dynamic spectrum v(t) derived for 


overdense loops (A5), which can be written as 


) (2IP) = y + -At 

e \ SV Q/ 2 A 


f o) j 


( 10 ) 


we perform a linear regression fit with the variables x h y h 
defined by 


Ti 


a 0 + (a! • x.) 
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Thus, for each time interval with a high-frequency cutoff v, •((,), 
we obtain from the linear regression coefficients ( a 0 , a x ) the two 
constants h 0 = a 0 X and v = a 1 X, which are the altitude h 0 = 
h(t = t 0 ) and the velocity v of an exciter that propagates in an 
overdense flux tube, specified by our density model (eq. [3]) 
and the parameters given in equation (A9). The so derived 
values for h 0 and v apply to fundamental (s = 1) plasma emis- 
sion, according to the arguments given in Appendix B. For 
harmonic (s = 2) plasma emission, the values for v and h 0 
would be a factor of s 2/p = 2 2/2,38 = 1.79 higher. 
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3.3. Selected Flare Events 

In the following we describe three selected flares, which show 
the existence of a high-frequency cutoff for three different types 
of radio emission: for a type IV event (Fig. 7), a decimetric 
pulsation event (Fig. 8), and dm type III bursts (Fig. 9). 

3.3.1. Type IV Event 

Figure 7 shows the radio data from Phoenix and SXR data 
from GOES. Since there are no HXR data available for this 
flare, we show the time derivative of the SXR data as a proxy 
for the HXR data (see discussion of the Neupert effect in 
Dennis & Zarro 1993). The first significant increase in the time 
derivative of the SXR flux (Fig. 7, bottom) occurs at 1200:50 
UT. The start of the GOES M3.6 flare was reported at 1 159 UT 
(Solar Geophysical Data), in coincidence with an Ha flare in 
NOAA Active Region 6718 at heliographic position S22E34. 
Enhanced decimetric radio emissi m was detected by the burst 


trigger of Phoenix at 1201 :26 UT, which is the start of the 
digitally recorded radio portion shown in Figure 7. This radio 
event was classified as type IV burst with fast-drifting fine 
structures and parallel-drifting bands (“ zebra ” patterns) in 
Isliker & Benz (1994; see Fig. 26 therein). 

The radio data show a persistent high-frequency cutoff 
throughout the event, which drifts slowly from 3 GHz initially 
(36 s after flare start in SXR) to 2.5 GHz (60 s after flare start). 
During this time interval there occur four episodes, of 1.2-4.5 s 
duration, showing a distinct high-frequency cutoff which 
exhibits a slow drift rate between —36 and —90 MHz s _1 in 
the frequency range of 2. 5-3.0 GHz. From these drift rates we 
infer with equation 4 upward-directed velocities between 60 
and 153 km s -1 . If we interpret the radio emission in terms of 
fundamental plasma emission, the observed frequency range 
corresponds to electron densities of n e = 0.9-1.2 x 10 11 cm -3 . 
These values appear to be remarkably close to the typical 
values inferred from the blueshift of the Ca xix line at the onset 



1200:00 


1210:00 


1 220:00 


Time [UT] 


Fig. 7.— Radio flux F R {top) and time derivative dF R /dt ( middle ) of Phoenix data from the flare on 1991 July 10, 1201 UT, in the frequency range of 2-3 GHz. The 
gradient-filtered dynamic spectrum ( middle ) shows fine structure with fast drift bursts and slowly drifting parallel bands. In the four time segments (marked with bars), 
drift rates of -90, -36, - 50, and -61 MHz s “ 1 are measured at the high-frequency cutoff - . The bottom frame shows the 0.5-4A SXR flux and time derivative from 
GOES. The flare start is indicated at 1200:50 UT. Note that the radio data extend over a interval of 25 s, while the GOES data are shown for a much larger interval of 
35 minutes. 
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of flares, with n e — 0.6-1. 0 x 10 11 cm 3 and v CE = 270 km s 1 
(Table 1). 

Our interpretation of this radio event in terms of chromo- 
spheric evaporation explains several observational properties: 
(1) the high-frequency cutoff” is consistent with free-free absorp- 
tion caused by the density enhancement (around 10 11 cm -3 ) of 
a chromospheric evaporation front, (2) the slow negative drift 
rate is consistent with an upward motion of the evaporation 
front with a speed of <150 km s _1 , and (3) the time delay of 
«36 s between the flare onset in SXR and the observed fre- 
quency cutoff at 3 GHz is consistent with the travel time of the 
chromospheric evaporation front upward to an altitude (»6.0 
Mm) corresponding to the plasma frequency level of 3 GHz in 
our model (see Fig. 2). This interpretation offers also a number 
of additional new insights: (1) the emission mechanism is 
uniquely identified as (fundamental) plasma emission, (2) the 
spectral peak of the radio emission is explained by the reduced 
free-free opacity due to the high temperature (J> 10 7 K) at the 


front of the upflowing plasma, (3) the parallel drifting “ zebra ” 
bands represent a wave pattern comoving with the density or 
temperature structure of the chromospheric evaporation front, 
perhaps associated with supersonic shock propagation, (4) the 
parallel displacement of the high-frequency cutoffs between the 
subsequent five episodes indicates either sequential injections 
or spatial separation of multiple chromospheric evaporation 
cycles. 

3.3.2. Decimetric Pulsation Event 

Radio and HXR data of the flare of 1992 October 5 are 
shown in Figure 8. The displayed time interval of 0923:23- 
0926:18 UT corresponds to the recording time of Phoenix after 
a burst trigger and is almost coincident with the time interval 
when BATSE/CGRO detected >25 keV HXRs, (i.e., 0923:25- 
0926:24 UT). SXRs detected by GOES start 43 s earlier than 
the HXRs, according to the first significant increase in the time 
derivative. An Ha subflare was registered in NOAA AR 7305 at 
heliographic position N12E41 (Solar Geophysical Data). 


v „ 

500 Fv 



0922:00 


0924:00 


0926:00 


Time [UT] 


Fig. 8. Flare of 1992 October 5, 0924 UT. The Phoenix radio data are shown background-subtracted (top), and gradient-filtered (middle). Both the low- and 

high-frequency cutoffs show a negative drift rate, - 14 and - 17 MHz s" 1 respectively, at the start of the fitted trajectory. The bottom panel shows HXR data > 25 
keV from BATSEjCGRO (histogram), SXR data at 0.5-4 A from GOES {monotonically increasing curve), and the time derivative of the GOES data (fluctuating 
curve). 
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The radio spectra show pulsating emission confined to the 
frequency range of 1. 0-2.0 GHz. Below 1 GHz the bandwidth 
of individual channels is 1 MHz (vs. 10 MHz above 1 GHz), to 
avoid interference. This enhances the receiver noise by a factor 
of ^/lO and reduces the sensitivity. Both the high- and low- 
frequency cutoff show a slow negative drift rate, which was 
determined to be — 14 and — 17 MHz s -1 at the start of the 
time segments indicated in Figure 8. This corresponds to velo- 
cities of 54 and 93 km s ~ 1 according to our density model. The 
evolution of the dynamic spectrum has a remarkable similarity 
with the simulated cases in Figure 5: (1) in the first 50 s after 
the flare onset the high-frequency cutoff of the radio emission is 
almost constant; (2) after 50 s (during the strongest impulsive 
rise of the HXRs) the radio start frequency shifts to a higher 
value and (3) then drifts slowly toward lower frequencies. The 
exciter velocity associated with the slowly drifting start fre- 
quency (93 km s -1 ) is lower than expected for a typical 
chromospheric evaporation front (300-400 km s -1 ), probably 
because the drift of the start frequency is evaluated over a time 
interval (104 s) that is much longer than elementary peaks in 
HXR (with timescales of 5-10 s) and therefore represents an 


envelope encompassing multiple injected evaporation fronts, 
similar to the case of multiple injections observed in Ca xix 
(Antonucci, Dodero, & Martin 1990). It is also interesting to 
note that the low- and high-frequency cutoffs of the radio emis- 
sion have an average ratio of 1:2, which corresponds to a 
density increase by a factor of 4. This is about the density 
increase that is expected for a chromospheric evaporation 
front, inferred from the densities of a flare loop at the start and 
end of the flare (see Table 1). 

3.3.3. Decimetric Type III Events 

In Figure 9 we show four short episodes of type III burst 
groups that exhibit a systematic, drifting high-frequency cutoff. 
These short episodes with durations of 10-15 s are more likely 
to display single evaporation cycles, while the other events 
shown in Figures 7 and 8 contain multiple evaporation cycles. 
The high-frequency cutoff of those longer events is likely to 
encompass a number of shorter episodes, which cannot always 
be disentangled properly. The drift rate of the high-frequency 
envelope is then systematically lower, which leads to an under- 
estimate of the velocity of the chromospheric evaporation front 
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Fig. 9. — Four episodes with type III emission and a drifting high-frequency cutoff from the flares on 1992 September 6 (top rows) and 1993 March 20 (bottom 
rows). For each case the background-subtracted (upper frame ) and the gradient-filtered (lower frame) radio flux is shown. The drift rates of the fitted curves at the 
high-frequency cutoffs are —82, —53, -56, and —42 MHzs' 1 . 
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TABLE 3 


Drift Rate and Velocity Measurements of High-Frequency Cutoffs 


Date 

^0 

(UT) 

D 

(s) 

v 0 

(MHz) 

dv/dt(t 0 ) 
(MHz s" 1 ) 

K 

(Mm) 

v 0 

(km s *) 

O 

I s 

i 

1991 Jul 10 

1201:28 

1.0 

2982 

-90.1 

6.0 

153 

11.0 

1991 Jul 10 

1201:31 

1.5 

3014 

-36.0 

6.0 

60 

11.2 

1991 Jul 10 

1201 :34 

3.3 

2962 

-49.9 

6.1 

85 

10.9 

1991 Jul 10 

1201:44 

4.5 

2765 

-61.3 

6.4 

120 

9.5 

1992 Sep 05 

1125:33 

13.0 

1878 

-15.3 

8.9 

60 

4.4 

1992 Sep 06 

0903:26 

49.0 

2079 

-14.0 

8.2 

46 

5.4 

1992 Sep 06 

1208:08 

12.0 

1786 

-82.4 

9.3 

360 

4.0 

1992 Sep 06 

1226:05 

11.0 

1484 

-53.0 

10.9 

325 

2.7 

1992 Sep 06 

1226:44 

13.0 

1157 

-5.4 

13.4 

53 

1.7 

1992 Oct 05 

0924:29 

104.0 

2106 

-16.7 

8.1 

54 

5.5 

1992 Oct 30 

1132:44 

40.0 

1327 

-7.7 

11.9 

58 

2.2 

1993 Mar 20 

1203:22 

15.0 

1607 

-56.7 

10.1 

300 

3.2 

1993 Mar 20 

1205:01 

23.0 

2292 

-115.3 

7.5 

318 

6.5 

1993 Mar 20 

1209:27 

45.0 

2205 

-62.0 

7.8 

184 

6.0 

1993 Mar 20 

1217:03 

42.5 

1233 

-5.7 

12.7 

49 

1.9 

1993 Mar 20 

1230:53 

10.0 

1298 

-41.7 

12.2 

327 

2.1 

1993 Mar 20 

1425:30 

25.0 

1704 

-27.4 

9.7 

130 

3.6 

1993 Oct 02 

0739:50 

30.0 

2147 

-23.1 

8.0 

72 

5.7 

1993 Oct 02 

0740:28 

15.0 

2104 

-70.6 

8.1 

228 

5.4 

1993 Oct 02 

0741 :26 

25.0 

1586 

-8.2 

10.3 

44 

3.1 

1993 Oct 07 

1228:52 

22.0 

1424 

-9.6 

11.2 

64 

2.5 

Average 


24.0 

1960 

-40.6 

9.2 

147 

5.2 



±23.3 

±590 

±31.7 

±2.3 

±114 

±3.1 

Average 

(5 < D < 15 s) 

12.7 

1616 

-46.4 

10.4 

236 

3.4 



±1.9 

±330 

±28.0 

±1.9 

±129 

±1.3 


(eq. [4]). Thus, the drift rates of the start frequency measured in 
shorter intervals such as those shown in Figure 9 are more 
likely to reveal the true velocity of resolved evaporation cycles. 
In the four cases in Figure 9 we measure drift rates of —82, 
— 53, —57, and —42 MHz s -1 . This corresponds to velocities 
of 360, 325, 300, and 327 km s - \ which are fully consistent 
with the velocity range of 300-400 km s -1 deduced from 
Ca xix lines. 

3.4. Statistics of Dynamical Parameters 
In Table 3 we summarize the measurements of all 21 
analyzed events possessing a high-frequency cutoff, with the 
parameters measured at the beginning (t 0 ) of the time intervals: 
the duration (£)), the frequency (v 0 ), and the drift rate ( dv/dt ). 
The density (n e0 ) is derived for fundamental plasma emission. 
The velocity ( v 0 ) and altitude (h 0 ) are inferred from the density 
model for overdense, type Ill-emitting flux tubes described 
above. The average velocity of all 21 measurements is v 0 = 147 
+ 114 km s“\ and the average height is h 0 = 9.2 ± 2.3 Mm. 
However, the longer episodes probably contain a systematic 
error due to the cumulative effect of multiple evaporation 


cycles. On the other hand, events which are too short (<5 s) 
have a larger uncertainty in the case of bursty emission. Thus, 
we consider the drift rates of start frequencies evaluated in 
intermediate time intervals, say in the range of 5-15 s, as most 
reliable. In Table 3 we find seven events with such durations, 
which imply a systematically higher value for the average 
upflow velocity; i.e., v CE = 236 ±129 km s -1 . 

3.5. Statistics of Timing at Flare Onset 
We examine the timing of the first detection of a high- 
frequency cutoff in radio with respect to the flare start detected 
in either hard X-rays or soft X-rays. From the 21 events 
analyzed in Table 3 we find seven episodes that occur at a flare 
start. The start times of the radio events, t 0 , are listed in Table 
4. The start times of >25 keV HXR emission, r HXR , and the 
HXR peak count rate (F H xr) are taken from the BATSE/ 
CGRO catalog (compiled at SDAC/GSFC). The earliest start 
of SXR emission, r SXR , is evaluated from the first significant 
increase in the time derivative of the 3-25 keV (0.5-4 A) GOES 
data. We exclude the event of 1993 October 7 from the sta- 
tistics because the detected flare start seems to be severely 


TABLE 4 

Timing of X-Ray and Radio Signatures 


*0 ( SXR J HXR ^HXR ( 0 f HXR *0 *SXR ^ o/ V 0 

Date (UT) (UT) (s) (counts- 1 ) (s) (s) (s) 


1991 JullO 1201:28 1200:50 ... ... ... 38 39 

1992 Sep 05 1125:33 1122:30 1122:03 95362 210 183 148 

1992 Sep 06 0903:26 0900:25 0901:45 99873 101 181 178 

1992 Oct 05 0924:29 0922:42 0923:25 24079 64 107 150 

1993 Mar 20 1203:22 1202:05 1202:08 25968 74 77 34 

1993 Oct 02 0739:50 0739:20 0739:29 316913 21 30 111 

1993 Oct 07 1228:52 1228:45 1228:55 1686 -3 7 175 

Average ... ... ... ... 94 103 110 

±71 ±67 ±61 
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affected by instrumental sensitivity limitations, since it has the 
lowest HXR count rate, being one to two orders of magnitude 
weaker than the other events. 

From the other six events we find for the earliest detection of 
a radio high-frequency cutoff an average delay of t 0 — t HXR = 
94 + 71 s with respect to the HXR onset, and t 0 — t SXR — 103 
+ 67 s with respect to the SXR onset. We compare these 
delays with the expected propagation delay At CE of the 
chromospheric evaporation front on its path up to the coronal 
height corresponding to the location of the radio source, using 
the parameters h 0 and v 0 from Table 3 (i.e., Af CE — h 0 /v 0 ). We 
find an average expected delay of At CE = 110 ± 61 s, which 
agrees well with the delay t 0 — t SXR measured after the onset of 
SXRs (103 + 67 s). 

4. DISCUSSION 

We discuss three reasons that substantiate evidence for the 
radio detection of chromospheric evaporation at decimetric 
wavelengths: (1) the propagation velocity (§ 4.1), (2) the elec- 
tron density (§ 4.2), and (3) the timing with respect to HXRs 
and SXRs (§ 4.3). Accepting the consequences of the chromo- 
spheric evaporation scenario as modeled here, several models 
of decimetric radio bursts require a revised interpretation that 
includes the dynamics of high-temperature flare plasmas 
(§ 4.4). 


4.1. Velocity Measurements 

The idea that the drift of the start frequency of decimetric 
radio bursts could be attributed to the chromospheric evapo- 
ration process was first considered in Aschwanden et al. (1993), 
where a velocity of a 300 km s -1 was inferred for the exciter 
associated with the high-frequency cutoff of a type IV event 
with quasi-periodic fine structure. 

In this study we investigated a variety of decimetric radio 
bursts (type IV, type Illdm, broadband pulsations), which all 
exhibit an upward drifting start frequency. We also found a few 
rare cases with downward drifting high-frequency cutoffs 
during short intervals (;$5 s) in some flares. Because HXR 
usually show major fluctuations on timescales of ^5 s 
(DeJager & DeJonge 1978), indicative of multiple energy re- 
leases in potentially different locations, it is likely that multiple 
chromospheric evaporation cycles occur. Such evaporation 
cycles have been observed on timescales of « 1 minute by 
Antonucci et al. (1990). For time intervals > 1 minute, it is 
likely that the averaged drift rate of the start frequency encom- 
passes multiple injections and thus yields only a lower limit to 
the speed of an individual evaporation front. We consider the 
measurement of the start frequency over time intervals of 5-15 
s more suitable for measuring the true propagation speeds of 
individual chromospheric evaporation cycles. Based on the 
statistics presented in Table 3, we find an average velocity of 
t> CE = 236 + 129 km for the shorter time intervals. Also, we 
have to be aware that the calibration of the velocity measure- 
ment used here is of statistical nature, relying on our model 
based on statistical values of type III drift rates. 

The drift rates of the detected high-frequency cutoffs are 
shown as a function of the frequency in Figure 10 for the 21 
analyzed events. The cutoffs fall into a band between 50 and 
350 km s _1 . We assumed a constant velocity in our fitting 
procedure, which is justified considering the small altitude 
range over which a fit extends. However, our velocity is mea- 
sured at a height of ^ 10 Mm, while soft X-ray blueshifts in 
Ca xix are weighted by the emission measure along the line of 
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Fig. 10. — Drift rates of type III bursts and high-frequency cutoffs. The 
average drift rate ( diamonds ) and standard deviation ( error bar ) of type III 
bursts are taken from Aschwanden & Benz ( 1986) (AB86) and Stahli & Benz 
(1987) (BS87). The upper line, labeled with v/c = 0.14, corresponds to the 
average drift rate of type III bursts, specified by the relation dv/dt = —0.01 
v 1 ' 84 of Alvarez & Haddock (1973). The lower straight lines, marking a velocity 
range of 50-350 km s -1 , encompass the drift rates measured at the high- 
frequency cutoffs of 21 flare episodes, labeled with date and start time on the 
right. 


sight. Since the centroid of the Ca xix emission is probably 
formed below 10 Mm, certainly during the first 30 s of the flare, 
the SXR-derived velocities apply to lower heights than in 
radio. 


4.2. Electron Density Measurements 
Radio emission at the plasma frequency provides the most 
direct measurement of the local electron density of a plasma. If 
radio type III bursts are observed, the only ambiguity in the 
local density measurement of the type III source is whether the 
radio waves are emitted at the fundamental (s = 1) or harmo- 
nic (s = 2) plasma frequency. This ambiguity in the frequency 
(v = sv p oc s^/nj causes a factor of four uncertainty in the elec- 
tron density. If we are to use the high-frequency cutoff of 
plasma emission as a density diagnostic of the chromospheric 
evaporation front, we have to resolve this ambiguity. Fortu- 
nately, there are three strong arguments that the observed 
high-frequency cutoff applies to fundamental plasma emission: 
(1) our model of overdense, type Ill-emitting flux tubes, based 
on the statistically observed type III drift rate, was found to be 
consistent with an average exciter speed of v lu /c = 0.14 in the 
case of fundamental plasma emission, which is consistent with 
the exciter speed v ni /c = 0.15 + 0.04 derived from absolute 
velocity measurements, while harmonic emission yields a less 
satisfactory agreement; (2) the mean velocity of chromospheric 
evaporation fronts inferred from our density model agrees 
better with the measurements of the Ca xix line in the case of 
fundamental emission (within a factor of 1.1) than for the case 
of harmonic emission (within a factor of 1.8); (3) the free-free 
opacity is able to produce a steep high-frequency cutoff for 
fundamental plasma emission in the 1-3 GHz frequency range, 
while no reasonable parameters can be found for harmonic 
plasma emission to produce the same effect (see simulated 
spectra in parameter space in Fig. 6). Thus, we conclude that 
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the slowly drifting high-frequency cutoff observed here is likely 
to be produced at the fundamental plasma frequency level. 
This gives us a very accurate tool to probe the electron density 
at the chromospheric evaporation front. For microwave type 
III bursts observed at 3.5 GHz (Stahli & Benz 1987) and at 
6.2-8.4 GHz (Benz et al. 1992), harmonic plasma emission was 
proposed based on the scale height argument derived from the 
observed drift rate, leading to a relatively low coronal tem- 
perature (0.5-0.7 MK) if the flux tube is assumed to be in 
isothermal equilibrium. In the present study we are forced to 
abandon the assumption of isothermal equilibrium in flaring 
flux tubes, especially in the environment of chromospheric 
evaporation. The high temperature («17 MK) of the hot 
upflowing plasma reduces the free-free opacity drastically 
(because of the T _3/2 dependence) and makes overdense flux 
tubes transparent to fundamental plasma emission, even at 
microwave frequencies. 

Our measurements indicate an electron density of n e = 
5.2(±3.1) x lO 10 cm -3 at the time of first detection of the 
chromospheric evaporation in radio, at an average height of 
h = 9.2 ± 2.3 Mm. From the emission measure of Ca xix, 
similar electron densities of 5.6-9.6 x 10 10 cm“ 3 (Table 1) are 
derived at the onset of the flare and typically increase by a 
factor of 4 during the flare. We have to take into account that 
the densities derived from the Ca xix line represent averages 
over the column depth below the chromospheric front, while 
the radio method measures the density directly at the top of the 
chromospheric evaporation front and thus returns a lower 
value. 

4.3. Timing of X-Ray and Radio Signatures 

The earliest detection of chromospheric evaporation indi- 
cated by a blueshifted Ca xix line generally coincides with the 
start of >25 keV HXRs (see Fig. 5 in Antonucci et al. 1993). If 
we take the start of >25 keV HXRs as a proxy for the first 
start of chromospheric evaporation, we can test the timing of 
the earliest detection in radio by the start time of the high- 
frequency cutoff. According to the statistics presented in Table 
4, we find the first radio detection of a decimetric high- 
frequency cutoff 94 + 7 1 s after the start of enhanced > 25 keV 
HXR emission or 103 + 67 s after the first significant increase 
in the rate of 3-15 keV SXR emission. The time difference 
between HXR emission and the time derivative of SXR emis- 
sion should be small for those flares where the Neupert effect 
(Dennis & Zarro 1993) applies and is indeed negligible here 
(9 s) compared with the delay times (94 or 103 s) of the radio 
onset. The onset of HXR or SXR corresponds to the launch of 
a chromospheric evaporation front down in the chromosphere, 
while the radio detection takes place at some height in the 
corona, and thus, a natural propagation delay is expected. 
Within the framework of our density model, we determined the 
altitudes h 0 and the velocities v 0 of the upflowing plasma self- 
consistently from the frequency and drift rate of the high- 
frequency cutoff. An interesting test is now to compare the 
extrapolated propagation times h 0 /v 0 with the actual delays. 
From the statistics presented in Table 4, we find an average 
extrapolated propagation delay oih 0 /v 0 — 110 + 61 s, which is 
in good agreement with the actual detection delay of t 0 
— t SXR = 103 ± 67 in radio with respect to the start in SXR. 
Thus, this temporal coincidence strongly supports a causal 
connection between the launch of a chromospheric evapo- 
ration front detected in SXR and the detection of a high- 
frequency cutoff in radio. 


Antonucci et al. (1990) identified successive injections of 
evaporating plasma on timescales of minutes, where each 
evaporation cycle showed initially a high velocity and 
decreased progressively during a cycle. Similarly, we identified 
multiple high-frequency cutoffs during a flare. The 21 episodes 
listed in Table 3 show a case with up to five episodes, which 
indicates multiple evaporation cycles. The intervals between 
successive episodes were found to be as short as 3 s (e.g., in 
1991 July 10 flare; Fig. 7), but it is not clear whether the 
multiple frequency cutoffs refer to temporally successive 
evaporation cycles or to different propagation delays resulting 
from different paths (e.g., adjacent flux tubes). 

4.4. Implications for Metric and Decimetric Radio Bursts 

Our chromospheric evaporation scenario predicts for times 
j>20 s after flare start (1) a slowly-drifting high-frequency (and 
possibly also a low-frequency) cutoff and (2) a significantly 
higher drift rate (than expected for the isothermal quiet 
corona) for decimetric bursts emitted at the fundamental 
plasma frequency. At first glance, a number of decimetric 
bursts can be found which match these morphological criteria: 
e.g., decimetric type III bursts (Figs. 3, 5, and 10 in Isliker & 
Benz 1994), type IV bursts (Figs. 24, 25, and 26 in Isliker & 
Benz 1994), and broadband pulsations (Figs. 28, 29, and 30 in 
Isliker & Benz 1994). Our chromospheric evaporation model 
sheds some new light on the interpretation of such burst types. 
As an example, we discuss the case of type III bursts with 
“ infinite drift rates.” 

The passage of a chromospheric evaporation front at fre- 
quencies of 200-500 MHz may have also been observed by 
Poquerusse (1994). In the examples of fast-drifting type III 
bursts shown in Figures 1 and 3 in Poquerusse (1994) there 
appears at «200 MHz a slowly upward drifting separatrix, 
which divides the normal-drifting metric type III bursts from 
the significantly faster drifting decimetric extensions. In our 
scenario this is naturally explained by the steeper density gra- 
dient at the chromospheric evaporation front which drifts 
slowly toward lower frequencies. Poquerusse (1994) proposes 
an alternative interpretation in terms of highly relativistic 
(v ul /c > 0.8) electron beams and suggests that these fast-drift 
bursts constitute a class separate from normal type III bursts. 
Since the chromospheric evaporation model provides a natural 
explanation for the transition of normal-drifting to fast-drifting 
type III bursts at decimetric frequencies, we do not see a strong 
reason to involve relativistic exciters. Poquerusse’s two classes 
of type III bursts may simply be renamed in type III bursts 
above and below the interface of a chromospheric evaporation 
front. 

5. SUMMARY AND CONCLUSIONS 
5.1. Observational Summary 

We analyzed the radio dynamic spectra in the 0. 1-3.0 GHz 
range from Phoenix, simultaneous HXR data from BATSE/ 
CGRO, and SXR data from GOES, from the years 1991-1993. 
We identified 21 time intervals during 11 different flares 
wherein the radio emission shows a high-frequency cutoff in 
the decimetric range, either in the form of a sharp spectral 
drop-off of continuum emission, or in form of a systematic 
change in the start frequency of bursty emission. The observa- 
tional results (see details in Tables 3 and 4) are as follows: 

1. The high-frequency cutoff during these 21 events shows 
a slow negative drift rate, with an average of ( dv/dt) 0 = — 41 
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± 32 MHz s \ measured at an average start frequency of 
v 0 = 1960 ± 590 MHz. 

2. The electron density of the average start frequency where 
the radio high-frequency cutoffs is first detected is n e0 = 
5.2( + 3.1) x 10 10 cm -3 (assuming fundamental plasma 
emission). This value is close to the value of n e = 5.6-9.6 
x 10 10 cm 3 inferred from the (blueshifted) Ca xix emission 
measure. 

3. The start time t 0 of the first detected high-frequency 
cutoff in a flare was found to be delayed by t 0 — £ H xr = 94 
+ 71 s with respect to the start of >25 keV HXR emission, 
and by t 0 — f SXR = 103 + 67 s with respect to the start of the 
3-25 keV SXR emission. 

5.2. Summary of Model-dependent Data Analysis 

We established a coronal density model for overdense, type 
Ill-emitting flare loops, consistent with the statistically 
observed drift rate of type bursts in the lower corona (power- 
law relation of Alvarez & Haddock 1973) and with an isother- 
mal, barometric model in the upper corona. This model was 
found to be consistent with a statistical exciter speed of v m /c = 

0.14 for type III bursts radiating at the fundamental plasma 
frequency, in agreement with the average value obtained from 
imaging observations, v lu /c = 0.15 + 0.04. Based on this model 
we find the following results: 

4. The exciter velocity associated with the slowly drifting 
high-frequency cutoffs corresponds to u 0 = 147 ± 114 km s -1 
for all episodes, or v 0 = 236 ± 129 km s -1 for episodes with 
5-15 s duration. In comparison, the velocity of the chromo- 
spheric evaporation front inferred from the blueshifted Ca xix 
line is 270 + 90 km s“ \ 

5. The altitude where the high-frequency cutoffs are first 
detected corresponds to h 0 = 9.2 ± 2.3 Mm. 

6. The extrapolated start time of an exciter starting in the 
chromosphere and propagating with velocity v 0 up to a 
coronal height h 0 is found to be h 0 /v 0 = 110 ± 61 s for the 
earliest episode in a flare. This value is close to the observed 
delay of the earliest detection of a radio high-frequency cutoff 
with respect to the onset of SXR, i.e., t 0 — t SXR = 103 ± 67 s. 

We interpret the slowly drifting high-frequency cutoff at 
decimetric wavelengths discovered here in terms of an opacity 
effect resulting from the upflowing hot flare plasma during the 
chromospheric evaporation process at the begin of the impul- 
sive flare phase. Evidence for this interpretation is given by the 
good statistical agreement between the parameters inferred 
from radio and Ca xix. 

5.3. Implications from Model Calculations 

We performed numerical calculations of the frequency- 
dependent free-free absorption for a propagating source radi- 
ating at the fundamental or harmonic plasma frequency. The 
configuration of the simulated model is an overdense flux tube 
embedded in an upflowing hot plasma, where the exciter of 
plasma emission is assumed to propagate inside the cylindrical 
flux tube. The flux tube has a diameter of 2-10 Mm and a 
coronal temperature T c = 1.5 MK, and the chromospheric 
evaporation front is simulated by a four-fold density increase 
and a temperature of T ce se 15 MK. The numerical results 
imply the following: 

1. The density gradient at the chromospheric evaporation 
front causes a local discontinuity in the plasma frequency that 


produces an almost infinite drift rate for type III— like bursts 
over the corresponding part of the radio spectrum (estimated 
to be within a factor of « 2 in frequency). 

2. The radio brightness is enhanced by a large factor of 
approximately (T ce /T c ) (3/2) ss 32 near the chromospheric 
evaporation front owing to the temperature dependence of the 
free-free opacity. 

3. Above the frequency v CE , the free-free radio spectrum 
rapidly falls off for fundamental plasma emission due to the 
density dependence of the free-free opacity. Harmonic plasma 
emission is not much affected in the same parameter range. 
Thus, the observation of a steep high-frequency cutoff implies 
fundamental plasma emission. 

4. The high-frequency cutoff for plasma emission is steeper 
for (1) lower chromospheric evaporation temperatures, (2) limb 
flares, or (3) small flux tube diameters. 

5. The high-frequency cutoff drifts toward lower fre- 
quencies, at a rate proportional to the velocity of the chromo- 
spheric front. 

Isothermal, hydrostatic equilibrium is an invalid assumption in 
the environment of chromospheric evaporation and should be 
abandoned when deriving the local electron density scale 
height from the drift rate of decimetric bursts, in particular for 
the fast-drifting decimetric bursts commonly detected during 
the impulsive phase of flares. 

5.4. Consequences for Future Work 

We have established a statistical correspondence between 
the SXR and decimetric signatures of chromospheric evapo- 
ration. The interpretation of decimetric radio emission during 
the impulsive flare phase should be reconsidered under this 
new aspect. The dynamics of the high-temperature plasma 
have been found to have a profound influence on the free-free 
opacity of decimetric plasma emission. In particular, we have 
demonstrated a new method for including the evolution of the 
upflowing hot flare plasma in physical models that can fit the 
temporal and spectral evolution of decimetric radio bursts. 
While this study is focused on statistical parameters to estab- 
lish the first evidence for a relationship between the pheno- 
mena observed in SXR and radio, we anticipate that future 
studies with simultaneous measurements in both wavelength 
regimes will allow a detailed modeling of individual flares. The 
BCS spectra from Yohkoh and simultaneous broadband 
dynamic spectra from Phoenix complement each other in an 
ideal way for this purpose. 

Studying chromospheric evaporation at radio wavelengths 
offers several advantages: (1) chromospheric evaporation can 
be detected from ground-based experiments, in contrast to soft 
X-ray measurements which require a spaceborne detector; (2) 
the plasma frequency measured in radio offers a very accurate 
density diagnostic of the upflowing plasma and complements 
the density measurement in SXR lines based on the line-of- 
sight integrated emission measure of the density-sensitive emis- 
sion lines; (3) the spectral shape at high-frequency cutoffs offers 
a sensitive temperature diagnostic: free-free absorption is espe- 
cially sensitive to admixtures of cooler plasma (T < 1 MK), 
while current SXR detectors are generally not sensitive in this 
temperature regime ; (4) no center-to-limb effect is expected for 
the velocity measurements inferred from drift rates at radio 
wavelengths, thus allowing velocity measurements even near 
the limb. This is not possible with the Doppler-shift method in 
Ca xix. On the other hand, blueshift measurements with the 



No. 2, 1995 


CHROMOSPHERIC EVAPORATION AND RADIO EMISSION 


1011 


Ca xix line offer a more accurate velocity measurement that 
can be used as a calibration of the density model. 

We found that the consequences of chromospheric evapo- 
ration are most effective for the decimetric range ( > 300 MHz). 
The expected effects at metric wavelengths are retarded by 
several minutes, the exciter velocity is slowed down, and the 
radiative signatures are perhaps diluted due to spatial disper- 
sion and merging of subsequent cycles. Possible related metric 
phenomena are type II bursts (shock waves), moving type IV 
bursts, or type I continua depressions (Bohme & Kruger 1982; 
Aurass et al. 1993). 

Finally, we would also like to point out that the radio detec- 
tion of upflowing, hot chromospheric plasma might be possible 
for stellar flares. Since coherent radio emission requires many 
fewer energetic particles than does incoherent (thermal) emis- 
sion at SXR wavelengths, the slowly drifting high-frequency 


cutoff of plasma emission should be easier to detect than a 
blueshift in SXR lines. Thus, the method described offers the 
possibility of detecting hot plasma upflows in stellar flares and 
to probe density models of stellar coronae. 
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APPENDIX A 

DENSITY MODEL FOR TYPE III-EMITTING FLARE LOOPS 

We construct a density model for flare loops which is (1) consistent with the statistically measured drift for type III bursts (eq. [2]; 
Alvarez & Haddock 1973) in the lower corona and (2) has a continuous transition to the isothermal, barometric atmosphere in the 
upper corona. This can be achieved by assuming a power-law function in the lower corona and an exponential form in the upper 
corona, i.e., 

h\~ p 


n e (h ) = i 


1 V 1 ! 


h < h x , 


n Q exp ( - - J h > h 1 , 


(Al) 


with a transition at height h x . The barometric model is constrained by the electron density n Q at the base of the quiet corona and the 
density scale height X. A smooth transition between the two regimes is obtained by requiring that the function and its first derivative 
be continuous at the transition height h = h v These continuity conditions determine the transition height h x and density n x = 
n e (h = h x ) at this height, i.e., 

h x = pX , n x = n Q exp ( — p) . (A2) 

For the lower corona (h < h x ) we calculate the dynamic spectrum v(t ) for the case of plasma emission. The observed frequency v is 
assumed to be emitted at the sth harmonics of the plasma frequency v p , and we can relate the frequencies to the electron densities 
according to (in cgs units) 


V = sv„ 


v„ = 


MY /2 - 8980 n l J 2 , v Q = MY = 8980 n^ 1 
\nm e J u \nm e J “ 


(A3) 


For an upward-moving exciter with constant velocity v, starting at time t 0 and height h 0 , the altitude h(t) can be parameterized as a 
function of time t by 

h{t) = h 0 + v(t - t 0 ) . 


Thus, for h < h x , 


v(t) = sv Q 


e h 0 + v(t — t, 
IP X 


,] 


p/2 


The time derivative of this dynamic spectrum yields the drift rate 

dv ev(sv Q )~ 2/p 


,d + 2 ip) 


dt 


2X 


(A4) 


(A5) 


(A6) 


Comparison with the empirical drift rate law of Alvarez & Haddock (1973), observed in the frequency range of 75-550 MHz, i.e., 
dv/dt = —Av a , with a = 1.84, and A MHz = 0.01 if the frequency is measured in MHz, or A = 9 x 10 -8 if measured in Hz, we find the 
following relations for the constants A and a : 


a = { 1 + - 
V 


A = 


ev{sv Q ) 


-2/p 


2X 


(A7) 
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This determines the remaining constants p and n Q in equation (Al): 

’-rri- mo) ~ 2 ■ (A8) 

Using the Baumbach- Allen model for the quiet corona, we find n Q = n^ A (h = 0) = 4.6 x 10 8 cm -3 . The scale height of the baro- 
metric model is determined by the temperature of the quiet corona ( T Q = 1.5 MK), given by the relation X — 4.6 x 10 9 T Q (MK). The 
numerical values for the constants in the density model n e {h) (eq. [Al]) are then 

p = 2.38 , X — 6.9 x 10 9 cm , h t = 1.6 x 10 10 cm , n Q = 4.6 x 10 8 cm -3 . (A9) 

Since this density model entirely relies on the statistical relation found by Alvarez & Haddock (1973), we check the consistency with 
recent measurements at higher frequencies. A comparison is shown in Figure 10. The centroid and width of drift rate distributions 
from 461 type III drift rates at 237 MHz, and from 294 type Illdm drift rates at 650 MHz (Aschwanden & Benz 1986), are found to 
be very close to the relation of Alvarez & Haddock (1973) (see symbols AB86 in Fig. 10). The average drift rate of 35 microwave type 
III bursts at 3.5 GHz was found to be 8100 ± 4300 MHz s _1 in Stahli & Benz (1987). This value (indicated with symbols BS87 in 
Fig. 10) is a factor of 4 below the expected value from the Alvarez & Haddock (1973) relation. 


APPENDIX B 

THE EXCITER SPEED OF TYPE III BURSTS 

The mean type III exciter velocity that is most consistent with our model inferred from the drift rate relation of Alvarez & 
Haddock (1973) can be obtained from equation (A8): 


j a ; 

[n Q (s 8980) 2 ] (1/p) . (Bl) 

e 

For the constants used above (eq. [A9]) we find u nI /c = 0.14 for fundamental (s = 1) plasma emission, or v m /c = 0.25 for harmonic 
(s = 2) plasma emission. In the following, we test which value is more consistent with measured exciter speeds of type III bursts. 

The values of measured exciter speeds of type III bursts compiled in Table 2 have been inferred with different methods that all 
have different systematic errors. Generally, there is a trend : a higher value was found in earlier studies using interferometric position 
measurements of type III sources, compared with interplanetary position measurements by later spacecrafts. Dulk et al. (1987) 
attribute this trend to the fact that (1) the apparent source location was not corrected for wave ducting, and (2) the density along the 
type III paths was overestimated in earlier measurements. It is beyond the scope of this work to evaluate the individual systematic 
errors of each measurement. Excluding the early measurements before 1970 which are believed to be systematically too high, we find 
an average value of y in /c = 0.15 + 0.04, by weighting the values according to the number of independently analyzed bursts. This 
mean value is remarkably close to the value found in our model for fundamental plasma emission, o in /c = 0.14. This constitutes also 
a new argument that most of these type III bursts are emitted at the fundamental plasma frequency. 
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